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Abstract 

QCD one-loop corrections to the semileptonic process e'^e" — ^"v^ u d 
are computed. We compare the exact calculation with a "naive" approach 
to strong radiative corrections which has been widely used in the literature 
and discuss the phenomenological relevance of QCD contributions for LEP2 
and NLC physics. 
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Introduction 



The measurement of VT-pair production at LEP2 will provide two additional pieces of 
information to our understanding of the Standard Model (SM) First, it will improve 
the determination of the W mass; second, it will probe the structure of triple gauge- 
boson couplings (TGC's). The mass of the W boson in the SM is tightly constrained. In 
fact an indirect determination of Mw can be obtained from a global fit of all electroweak 
data. The fit gives 

Mw = 80.359 ± 0.05ll|^:|J^^GeV (1) 

where the central value correspond to Mh = 300 GeV and the second error reflects 
the change of M-^ when the Higgs mass is varied between 60 and 1000 GeV. A more 
precise determination of Mw will provide a stringent test of the SM. A disagreement 
between the value of Mw derived from the global fit and the value extracted from direct 
measurement would represent a major failure of the SM. Alternatively, an improvement 
in the value of the W mass can significantly tighten present bounds on the Higgs mass. 

Two methods have been singled out as the most promising [0]. The first one is 
based on the rapid increase of the total cross section at threshold. The second method 
relies on the direct reconstruction of the mass from the hadronic decay products of the 
W using the decay channels 

W^W- qqiu (2) 

W+W- ^ qiqiq2q2 (3) 

where i = e, fj,. 

Preliminary studies indicate that the direct measurement will provide a more precise 
determination of the W mass than the threshold method for which a smaller number 
of events will be available. Combining all decay channels, an accuracy of about 35 
MeV is expected from direct reconstruction, while the ultimate precision attainable at 
threshold is estimated to be about 100 MeV. 

In WW production triple gauge-boson couplings, which play a central role in non- 
abelian gauge theories, appear already at tree level and can be studied in much more 
detail than at lower energies [0]. At LEP2 TGC's will be mainly probed using angular 
distributions of the W^s and of their decay products. At the higher energies available 
at the Next Linear Collider (NLC) it will be possible in addition to study TGC's using 
the energy dependence of the total production cross section, since non-standard TGC's 
in general lead to a cross section which increases with energy and therefore eventually 
violates unitarity. 

On-shell VT-pair production is described at tree level by a set of three diagrams, 
labeled (e) and (f) in fig. 1. At the lowest level of sophistication, one can attach 
to them the on-shell decay of the W^s. In this case all one-loop electroweak and 
strong radiative corrections are known 0, |^, This approach is gauge invariant, but 
all effects of the finite width of the W and all correlations between the two decays 
are neglected. In order to take these features into account one can consider the set 
of diagrams mentioned above with off-shell W^s, which decay to four fermions. In 
the literature the corresponding amplitude is called CC3. CC3 however is not gauge 
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invariant. A gauge independent description of WW production requires in the unitary- 
gauge ten (twenty) diagrams for the semileptonic channel (|^) when £ = /i (e) and eleven 
diagrams for the hadronic channel (^. These amplitudes are known as CCIO (CC20) 
and ceil respectively. Numerically the total cross sections obtained from CCS and 
CCIO or ceil at LEP2 and NLC energies differ by a few per mil. Larger discrepancies, 
of the order of several per cent, are found when comparing CCS with CC20. 

In order to extract the desired information from WW production data, theoretical 
prediction with uncertainties smaller than those which are foreseen in the experiments 
are necessary. In particular this means that radiative corrections have to be under 
control. 

In this letter we present the complete calculation of QCD corrections to CCIO. In 
most studies they have been included "naively" with the substitution Tyy ^^{1 + 
2as/7r/3) and multiplying the hadronic branching ratio by (l + a^/vr). This prescription 
is exact for CCS when fully inclusive quantities are computed. However it can only be 
taken as an order of magnitude estimate even for CCS in the presence of cuts on the 
jet directions, as discussed for the hadronic channel in |^ [|. It is well known that 
differential distributions can be more sensitive to higher order corrections than total 
cross-sections in which virtual and real contributions tend to cancel to a large degree. It 
is therefore necessary to include higher order QCD effects into the predictions for WW 
production and decay in a way which allows to impose realistic cuts on the structure 
of the observed events. 




Figure 1: Tree level diagrams for e~^e — > /i u d. The dashed lines are Ws. 



^ The impact of QCD corrections on the angular distribution of the decay products oi a W and 
their application to on-shell VF-pair production is discussed in ref. M 
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Figure 2: Basic combinations of loop diagrams. All virtual QCD corrections to electroweak 
four fermion processes can be computed using these three sets. The quark wave functions 
corrections are not included because they vanish, in the massless limit, using dimensional 

regularization. 



Calculation 

One-loop virtual QCD corrections to e+e^ ^^v^ u d are obtained by dressing all 
diagrams in fig. 1 with gluon loops. Defining suitable combinations of diagrams as in 
fig. 2 one can organize all contributions in a very modular way, as shown in fig. 3. 
Note that the first and the fourth contributions in fig. 3 can be obtained from each 
other by multiplying by -1 and interchanging momenta and helicities of u and d. 

The calculation has been performed using standard Passarino-Veltman techniques 
H and dimensional regularization for ultraviolet, coUinear and soft divergences. With 



the help of the Symbolic Manipulation program FORM |]I0[ , all tensorial integrals have 
been reduced to linear combinations of scalar loop functions. One needs to calculate 
one four-point, one two-point and three three-point basic functions with different input 
momenta, so that twenty independent scalar loop functions contribute to the cross 
section. 

Having classified loop corrections as in fig. 2, one easily convinces oneself that ex- 
actly the same ingredients appear in the computation of 0{as) virtual corrections to 
any electroweak four-fermion process. In fact, due to the color structure, gluons con- 
necting different spinor lines in the final state start contributing at 0{a1). Therefore, 
all that is required in order to extend our results to the calculation of QCD corrections 
for all possible electroweak four-fermion final states is the computation of the real gluon 
emission amplitudes. We plan to pursue this program in the near future. 
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Figure 3: One loop gluonic corrections to e~^e ^ u d. 

The real emission contribution for e^e~ ^'v^ u d can be obtained attaching a 
gluon to the quark hne of the diagrams shown in fig. 1 in all possible positions. This 
results in twenty-four diagrams. The required matrix elements have been computed 
using the formalism presented in ref. with the help of a set of routines (PHACT) 
1131 which generate the building blocks of the helicity amplitudes semi-automatically. 

If we write the full NLO cross section in the form 

ONLO = J^da"^ + J^da^ , (4) 

where we assume that all ultraviolet divergencies have been canceled by renormaliza- 
tion, the real (R) and virtual (V) contributions are still separately singular in four 
dimensions because of soft and collinear singularities, while the sum is finite. 

In order to be able to integrate separately the real and virtual part one has to 
explicitly cancel all singular contributions in each term in a consistent way. To this 
aim we have used the subtraction method: the full expression is rewritten in the 
form 

aj^LO = ^ [da^ - da^] + J/^^ + ^ ^a^ ■ (5) 

The subtraction term da^ must have the same pointwise singular behaviour as the 
exact real emission matrix element in order to cancel soft and collinear divergencies. 
It must also be possible to integrate da^ analytically in d dimensions over the one 
parton subspace which generates the singularities. The result of this integration is then 
summed to the virtual contribution producing a finite remainder that can be treated in 
four dimensions. Benefits of this method are twofold. First, an exact result is obtained 
and no approximation needs to be taken; second, all singular terms are canceled under 
the integration sign and not at the end of the calculation, leading to better numerical 
accuracies. This is especially relevant for the present case, since we are aiming for high 
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parameter 


value 


Mz 


91.1888 GeV 


Tz 


2.4974 GeV 


Mw 


80.23 GeV 


Tw 


3G'^M3,/(v/87r) 


a-^ = a-\2Mw) 


128.07 


Gp 


1.16639- 10^^ Gev-2 


sin^ Ow 




as 


.111 (.123) 


VcKM 


1 



Table I: Input parameters. 



precision results, with errors of the order of a per mil. We have found it particularly 
convenient to implement the recently proposed dipole formula |]ll| . These are a set of 
completely general factorization expressions which interpolate smoothly between the 
soft eikonal factors and the collinear Altarelli-Parisi kernels in a Lorentz covariant way, 
hence avoiding any problem of double counting in the region in which partons are both 
soft and collinear. 



All integrations have been carried out using the Montecarlo routine VEGAS ||T4|. 

An important ingredient for accurate predictions of ly-pair production is the effect 
of electromagnetic radiation. In the absence of a calculation of all 0{a) corrections to 
four-fermion processes, these effects can only be included partially. In contrast with 
LEPl physics a gauge invariant separation of initial and final state radiation is not pos- 
sible. On the other hand, the leading logarithmic part of initial state radiation is gauge 
invariant and can be included using structure functions. The non-logarithmic terms 
however are unknown. In order to assess the influence of QCD corrections, we are inter- 
ested in a comparison with the results obtained in the Workshop on Physics at LEP2. 
Hence we have decided to employ as much as possible the parameters adopted in the 
"tuned comparisons" |]T5[ which provide the most extensive collection of results. There- 
fore we have used the (3 prescription in the structure functions, where (3 = ln(s/m^) — 1. 
In the same spirit we have not included Coulomb corrections to CC3, which are known 
to have a sizable effect, particularly at threshold. They could however be introduced 
with minimal effort. 
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Born 


Exact QCD 


"naive" QCD 


161 GeV 


.24962 ± .00002 


.24760 ± .00002 


.24790 ± .00002 


175 GeV 


.96006 ± .00007 


.94519 ± .00007 


.94613 ± .00007 


190 GeV 


1.184003 ± .00009 


1.16681 ± .00009 


1.16766 ± .00008 


500 GeV 


.46970 ± .00006 


.47109 ± .00007 


.46131 ± .00006 



Table II: Cross sections in pb with canonical cuts (see text). 



Results 

In this section we present a number of cross sections and of distributions for e"'"e~ 
H'u^ qi 7j2- In all cases we sum over the two possibilities (gi, ^2) = {u, d) and (gi, ^2) = 
(c, s). The input parameters used in our calculation are given in table I. At LEP2 
energies we have used as = .117 as in while at the NLC we have adopted as = .123 
in order to conform to the choice made for the Joint ECFA/DESY Study: Physics and 
Detectors for a Linear Collider. Initial state radiation is included in all results. 

Previous studies [|15| have shown that the differences between the total cross sections 
obtained from CCIO and those obtained with CC3 are at the per mil level. Much larger 
effects have been found in observables like the average shift of the mass reconstructed 
from the decay products from the true W mass. If S-_ and s+ are the invariant masses 
of the /i~z/^ pair and of the hadronic system, respectively, the standard definition is: 



(AM) = ± ^'""C " da. (6) 




This quantity vanishes in the zero width approximation and provides a useful estimate 
of the influence of various physical processes on the relationship between the measured 
value of the W mass and its actual value. 

For LEP2 we have adopted the so called ADLO/TH set of cuts: 

• > 1 GeV; 10° < 9^ < 170° 

• < 5° 

Furthermore: 

• the energy of a jet must be greater than 3 GeV; 

• two jets are resolved if their invariant mass is larger than 5 GeV; 

• jets can be detected in the whole solid angle. 
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For the NLC we have adopted the NLC/TH set of cuts which differ from the 
ADLO/TH set in that a minimum angle of 5° is required between a jet and either 
beam and that two jets are resolved if their invariant mass is larger than 10 GeV. Both 
set of cuts will also be referred to as "canonical" in the following. 

The assumption that each final state particle corresponds to a jet must be abandoned 
when going from LO to NLO calculation. Starting from the final state partons, it 
is necessary to define jets using an infrared safe procedure. Only in this case the 
cancellation of infrared and collinear singularities between virtual and real corrections 
can take place and meaningful results can be obtained. In the present case it is natural 
to define jets using the ADLO(NLC)/TH cuts. Following this prescriptions we have 
merged into one jet those parton pairs whose invariant mass was smaller than 5(10) 
GeV. Furthermore partons with energy below 3 GeV have been merged using the JADE 
algorithm. Having identified jets, we have checked whether they passed the canonical 
cuts. All events with two or three observed jets have been retained in our plots and 
cross sections. 

In fig. 4 we present the normalized distribution of the angle between the muon and 
the closest jet at the three energies at which data will be taken at LEP2, ^/s = 161, 175 
and 190 GeV. The corresponding cross sections are given in table II together with the 
results at Born level and those obtained with "naive" QCD (nQCD) corrections. The 
events in fig. 4 pass all ADLO/TH cut with the exception of the minimum Notice 
that imposing the latter cut corresponds to discarding only the leftmost bin in the 
plot. According to the ADLO/TH prescription, jets can be measured over the full solid 
angle. Therefore, apart from the consequences of the increase in the W width, the 
only possible effects of QCD radiation in the semileptonic channel are related to the 
increased probability of a jet to be close to the charged lepton. With increasing collider 
energy, the two VF's tend to fiy apart with larger relative momentum and therefore 
the probability of a jet to overlap with the lepton decreases. This behaviour is clearly 
visible in fig. 4 and is confirmed by fig. 6a which shows that at y/s = 500 GeV jets 
are typically well separated in angle from the charged lepton and that the differences 
between the exact distribution and the one with "naive" corrections is confined to very 
large angles. From fig. 4 and 6a it is also apparent that, as expected, exact QCD 
predicts smaller minimum angle between jets and the charged lepton. nQCD results 
for this observable become closer to the exact NLO distribution as the energy increases. 

Table II shows that QCD corrected cross sections for e+e~ fJ'^^ii Qi Q2 are between 
1%, at ^ = 161 GeV, to 2% smaller, at = 190 GeV, than Born prediction at LEP2 
with ADLO/TH cuts. In the absence of cuts the exact NLO result for the total cross 
section and the "naive" implementation of QCD corrections agree to better than 10~^ 
and are indistinguishable within statistical integration errors for all energies studied 
in table II. This is expected since the two inclusive cross sections differ by corrections 
of order ag/n to terms of relative order 10^^. After cuts the two sets of results at 
LEP2 differ by about 10~^, while statistical errors are smaller than 10~^. Therefore, 
not surprisingly, there is clear evidence that nQCD fails even for cross sections when 
phase space for final particles is limited by cuts. However, with LEP2 canonical cuts 
the difference is much smaller than the projected experimental accuracy. Because 
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of the additional cut on the minimum angle between jets and initial state leptons, 
measurements at the NLC are far less inclusive than at LEP2. The effect is further 
amplified by the high-energy behaviour of ly-pairs which tend to be produced at 
smaller angles that near threshold. Therefore at ^/s = 500 GeV, with NLC/TH cuts, 
the exact NLO cross section is about 2% larger than what is obtained with nQCD, 
a difference which is much larger than the expected experimental uncertainties. It is 
amusing to note that in this case the exact NLO result is very close, about .5% larger, 
to the tree level cross section than to the nQCD prediction, which is approximately 
1.5% smaller than the Born result. The decrease in cross section due to the larger 
width of the W is compensated by the extra radiation which leads to a larger number 
of events with two or more visible jets. 

The mass shift (AM) at y/s = 175 GeV, with canonical cuts, is found to be 

(AM)jvLo = -0.6383 • 10"^ ± 0.0002 • 10"^ (7) 

which happens to be in excellent agreement with the nQCD result {AM)nQCD = 
—0.6381 ■ 10~^ ± 0.0002 ■ 10~^. This is to be compared with the tree level result of 
{AM)Born = -0.6219 ■ 10-2 ± 0.0002 • lO'^. Therefore the large effect, about 2.5%, 
of QCD corrections on (AM), which was suggested by the "naive" approach |]T5[, is 
confirmed by our calculation and placed on a solid footing. 

Fig. 5 shows the distribution of the minimum angle between any jet and either beam 
at ^/s = 175 GeV. Even though the ADLO/TH set does not include a requirement on 
the 9j angle, such a cut is included in the experimental studies for the hadronic channel, 
and is part of the NLC/TH set. We have separated the contribution of two-jet and 
three-jet final states making it possible to estimate the effect of different angular cuts 
on the cross section. 

In fig. 6 two distribution at y/s = 500 GeV are presented. Fig. 6a shows the 
distribution of the angle between the muon and the closest jet, while fig. 6b shows the 
distribution of the minimum angle between any jet and either beam. Fig. 6a makes it 
clear that essentially all events would pass any reasonable isolation cut for the In 
fig. 6b we separate again the contribution of two-jet and three-jet final states so that 
the effects of non-canonical cuts can be judged from the plot. 

Conclusions 

We have described the complete calculation of QCD radiative corrections to to the 
semileptonic process e'^e~ —>■ fi~u^ u d which are essential in order to obtain theoreti- 
cal predictions for VT-pair production with per mil accuracy. The amplitudes we have 
derived are completely differential, and realistic cuts can be imposed on the parton level 
structure of the observed events. It has been shown that the "naive" implementation 
of QCD corrections fails in the presence of cuts on the direction of final state jets. 
The error is typically at the per mil level at LEP2. However at ^/s = 500 GeV, with 
NLC/TH cuts, the discrepancy is about 2%, much larger than the expected experimen- 
tal precision. QCD corrections substantially increase, by more than 2%, the average 
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shift between the W mass measured from the decay products and the actual value of 
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Fig. 4: Distribution of the angular separation of the /j, from the closest jet at \/s = 
161, 175, 190 GeV. All ADLO/TH cuts with the exception of that on e'^L]^^ arc applied. 
The continuous histogram is the exact NLO result while the dashed histogram refers 
to nQCD. 
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Fig. 5: Distribution of the minimum angular separation between any jet and either 
beam at ^/s — 175 GeV with canonical cuts. The continuous, dotted and dot-dashed 
histograms are exact NLO results while the dashed histogram refers to nQCD. 
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Fig. 6: Distribution of the angular separation of the /i" from the closest jet (a) and 
of the minimum angular separation between any jet and either beam (b) at y/s — 500 
GeV. The continuous, dotted and dot-dashed histograms are exact NLO results while 
the dashed histogram refers to nQCD. 



